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A linear external perturbation is introduced in the action of the partition function of 
the random matrix model of RNA [G. Vernizzi, H. Orland and A. Zee, Phys. Rev. Lett. 
94, 168103 (2005)]. It is seen that (i). the perturbation distinguishes between paired and 
unpaired bases in that there are structural changes, from unpaired and paired base structures 
(0 < a < 1) to completely paired base structures (a — 1), as the perturbation parameter a 
approaches 1 (a is the ratio of interaction strengths of original and perturbed terms in the 
action of the partition function), (ii). the genus distributions exhibit small differences for 
small even and odd lengths L, (iii). the partition function of the linear interacting matrix 
model is related via a scaling formula to the re-scaled partition function of the random matrix 
model of RNA, (iv). the free energy and specific heat are plotted as functions of L, a and 
temperature T and their first derivative with respect to a is plotted as a function of a. The 
free energy shows a phase transition at a = 1 for odd (both small and large) lengths and 
for even lengths the transition at a = 1 gets sharper and sharper as more pseudoknots are 
included (that is for large lengths). 

PACS numbers: 02.10.Yn, 87.14.gn, 87.10.-e, 02.70.Rr 



'Electronic address: ittygarg@physics.du.ac.in 
t Electronic address: ndeo@physic s.du.ac.in| 



2 



I. INTRODUCTION 

Ribonucleic acid (RNA) is the only known bio-molecule that plays the dual role of being a 
carrier of genetic information (from Deoxyribose- nucleic acid (DNA) to proteins) and an enzyme 
in important biological reactions The secondary structure of RNA, excluding pseudoknots, has 



been a central subject of study for understanding folding of RNA [2j, |3|. However, experimental 
studies [J] have shown that tertiary structures of RNA play pivotal role in determining their 
biological functions. Therefore, in order to understand the role that RNA's play in biological 
processes it becomes essential to understand their folded conformations. The RNA folding problem 
has been studied using a number of statistical models In this work we study the effect of 
introducing an external interaction in a random matrix model of RNA as RNA in a cell participates 
in important genetic processes like transcription and translation, is surrounded by ions, comes 
in contact with cell walls, other molecules and organelles. These models are inspired from a 
theoretical model of RNA which maps the RNA folding problem onto a large N matrix field theory 
[(J where N — > oo limit gives secondary structures and tertiary structures (pseudoknots [7]) are 
obtained by finding corrections to 1/N 2 terms in the partition function. The partition function 
can be expanded in powers of jhg which gives a topological expansion, originally observed by 
G 't Hooft [§] where planar and tertiary structures can be identified as terms with g = and 
g > 1 respectively (g is the genus of RNA structures). Thus random matrices give a natural way 
of arranging the RNA structures according to their topology which is crucial to study the genus 
distribution functions of the model. Here it is natural to introduce the chemical potential \i > 
via the activity k = exjr^ 1 = 1/N 2 , then k — > gives the planar diagrams and k — > oo yields 
the maximum possible average g. Thus the topological chemical potential controls the number 
of pseudoknots in folded conformations (Eq. (6) in [f| can be written as Eq. (4) in Q]), hence 
connecting N in these matrix models to a physical quantity fi. This physical connection has been 
elaborated on recently in [l(J 



The random matrix model of RNA in |ll|] enumerates all planar structures along with the 
pseudoknots by making the following simplifying assumptions in the model of [6] : (i) base pairings 
are not complementary i.e., they do not depend on the type and location of a monomer in the 



3l(]. This 



chain and have the same strength v (homopolymer) and (ii) the chain is infinitely flexible 
enables us to count the RNA structures with pseudoknots and also to observe how they evolve as 
the length of the chain is increased. Therefore, it addresses an important fundamental question 
in the prediction of RNA structures i.e., exact combinatorics of RNAs with pseudoknots 15] 32]. 
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However, the number of structures that have been discovered are a very small subset of the vast 
number of structures given by the matrix model. Therefore, in order to study the effect of the 
surrounding environment on RNA and its conformations, addition of a perturbation in the action 
of the random matrix model partition function of [111 ] and to observe the changes is important. 

As a first step, section 2 discusses an external linear interaction in the action of the parti- 
tion function of the random matrix model in 11)]. The corresponding distribution functions and 



thermodynamic properties are studied in sections 2 and 3. 

II. MATRIX MODEL OF RNA WITH A LINEAR EXTERNAL INTERACTION 

Here the matrix model of RNA with an external perturbation is presented and discussed. A 
perturbation term U(4>), which can be chosen to be a polynomial with linear, cubic, quartic and 
other higher order terms in <j), is added in the action of the partition function of the matrix model 
in [11 ] (keeping all the assumptions of the RNA matrix model in [ll| ) . The partition function of 
the matrix model with an external interaction is the integral 

Z L {N) = [ n^e-f ^rt^V^WilVna + W. (1) 

L( - ' J 1=1 1=1 

As a first step U{4>) is considered to be a linear perturbation of the form ^2i = i{Wi)~ 1 Tr<pi 
where W» acts on each </>j. So the partition function becomes 

Z L (N) = —L- [ n^e-f ^(v-Jo^^E^W)- 1 ^ 1 Tr ^ (1 + (2) 

L ^ ' i=l i=l 
The 4>i's are i = 1, L independent (N x N) random hermitian matrices. The matrix elements 
of (j>i in both the quartic and perturbation terms are taken from the same random Gaussian 
distribution with zero mean and unit variance. A simplification is made, W% = w, which makes the 
interaction act uniformly on each V is an (L x L) interaction matrix (jj containing interactions 
between different </>j's which is considered to be Vij = v = e~@ e (here = 1/kgT with kg as the 
Boltzmann constant and e as the base specific pairing energy which is assumed to be same for any 
base pairing i.e., a homopolymer which is also infinitely flexibility) and the observable Y\i(^ + fa) 
is an ordered matrix product over fa's which evaluates to a polynomial of order L. The observable 
ensures that the diagonal elements Vu of the interaction matrix V in the quartic term of the action 
do not appear in the partition function Zl{N). Al(N) is the normalization constant given by 
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A L {N) = / Yld^e-^^'^^e-^^- 17 ^. (3) 
J i=l 

The matrix model partition function Eq. (2) can be written by raising the term jjTr Yii=i 
in the exponent as exp log [jv r - rr n<=i( 1 +^»)] anc [ writing the partition function with a complicated 
action : (=f E£=i(^%-2><^ ~ # EtiCWi)- 1 :^ - log [^Tr nf =1 (l + &)])• In this rep- 
resentation, performing a perturbative expansion near a minimum of the action by keeping the 
quadratic part of the action in the exponential and expanding the non-quadratic terms, results in 
computing series of polynomial moments of Gaussian integrals which are represented diagrammat- 
ically by the Wick's theorem. The normalization constant Al(N) under simplified assumptions 
Vij = v and Wi = w can be written as A L (N) = exp^ Tr{ -^ /Il£=i d® i exp~^ Tr ^ i ( v ~^ ij ' s '^ where 
<E> is defined as $j = {(^i + VijW^ 1 ). Carrying out a series of Hubbard Stratonovich Transformations 
in the multi matrix integral of Eq. (2) and making a redefinition a' = (a + ^) = (a + a), where 
a = — is defined as the ratio of strength of interaction between vertices to the strength of the 

to ° ° 

linear external perturbation, reduces the integral to be over a single (N x N) matrix a 



where Ri, a {N) = f daexp~^ Tr ( a ^ . The integral in Eq. (4) is a Gaussian integral with an 

observable which gets shifted by an amount a as a result of the redefinition of a. For comparison, 
I I 

the model in [17[ is considered where potentials such as U(x) = g\x + ^x 2 + |x 4 are studied. The 
partition function for this model can be written by keeping (gix + ^x 2 ) in the exponent and pulling 

down expi Trx4 as a series, expi TrxA = (1 + |Trx 4 + ^(| ) 2 (Trx 4 ) 2 + ), as done here. Then 

(<7i x + t;X 2 ) can be redefined so that matrix x shifts to x' to get a Gaussian in the exponent. The 
series or observable as a function of x' is, exp J r x 2g i = [1 + |Tr(V — 2^") 4 + <jr(f ) 2 (Tr(x' — 

2^j-) 4 ) 2 + ]. This is analogous to the observable jjTr(l + a' — a) L in Eq. (4). 

The spectral density /OjVa(A) of a Gaussian matrix model defined at finite N is given by 



Making use of the identity /_ °° d\pN ja (\) = 1 in Eq. (4), the partition function can be written 

as 
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Z L ,a( N )= I d\p N>a (X){l + X-a) L . 

J —oo 

Define G(t,N,a) to be the exponential generating function of Z^^N) 



(6) 



G(t,N,a) = J2 Z LA N ) 



L=0 



TV 



G(t, N, a) can now be written using Eq. (6) as 



G(t,N,a) = / d\p Nja (\)J2 



t L (l + \-a) 



CO 



L=0 

d\p N)a (\)e t{1+x - a \ 



The form of spectral density pN,a (A) from Q,Q 



jV(A) 2 N-l IT I \ [~N \ 



(7) 



(8) 



(9) 



where H2k{^\J~^) represents Hermite polynomials, is used in Eq. (8) to write G(t,N,a) as 



G(t,N, a) 



(10) 



Completing the square in the above equation and using a standard result of integration over 
Hermite polynomials [20j, f*™ dx e~( x ~ y } 2 H n (x) = ^y n 2 n , Eq. (10) solves to 



G{t,N,a) = e w+*( 1 -«) 
Combining Eq. (7) and Eq. (11), G(t,N,a) can be written as 



N ^ ~ {k+1) k\N k 

k=0 



(11) 



G(t,N,a) = y j Z L AN)- = e^ 1 



t &+t(l-a) 



L=0 



fc=0 



(12) 



The parameter a in Eq. (12) appears only in the exponent. The partition functions for different 
L and a can be obtained by comparing powers of t on both the sides of Eq. (12) (Table I) and is 
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FIG. 1: The figure shows a 3D plot of the partition function Zl,c(N) of the linear interacting matrix model 
of RNA as a function of a(= 0,0.25,0.5,0.75,1) and L (upto 20) in (a). The plot in (b) magnifies the 
a = 0.75, 1, 1.25 region. At a = 1, Zl^N) has a finite non-zero value for even lengths whereas it is zero 
for odd lengths. 

shown in Fig. 1(a) for L < 20 and < a < 2 (raised in intervals of 0.25). For a = 1, only even 
L partition functions are non-zero whereas the odd L partition functions vanish (Fig. 1(b)). For 
completeness we give the explicit dependence of partition function Zl^{N) on N in terms of the 
topological parameter, genus g, as Zl j0! (N) = Y^=o a L,g,aj^ [1 11 ] - Here the coefficients aL,g,a give 
the weighted number of diagrams at a given L, genus g and a. Also, the total weighted number of 
diagrams are given by M a = Zl^{N = 1), for a particular L and a, independent of the genus. 

It is found on doing the matrix integrals, that the partition function and weight of any config- 
uration (Feynman diagram) are both positive for the interaction parameter a when < a < 1, 
listed in Table I. But when a > 1, the partition function and weights of Feynman diagrams be- 
come negative for odd lengths of the chain, however they still remain positive for the even lengths. 
Since each Feynman diagram represents a conformation of RNA which must have a non-negative 
weight, Feynman diagrams of the matrix model with linear external interaction correspond to RNA 
structures for < a < 1 for odd L and for all values of a for even L. 



The general form of the partition function for the matrix model with linear external interaction 
can be obtained from Eq. (2) by writing the partition function in terms of the variable $ (introduced 
in section 2) as 



A. General form of Z L ^ a (N) 
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TABLE I: The Table lists partition functions Zi,. a (N) for different lengths L at any a for the linear inter- 
acting matrix model of RNA. 

L Z L , a (N) 

1 (1-a) 

2 (l-a) 2 +v 

3 (l-a) 3 + 3v(l-a) 

4 (l-a) 4 + 6v(l-a) 2 + 2v 2 + v 2 /7V 2 

5 (1 - a) 5 + 10w(l - a) 3 + 10v 2 (l - a) + 5v 2 (l - a)/N 2 

6 (1 - a) 6 + 15w(l - a) 4 + 30v 2 (l - a) 2 + 15v 2 (l - a) 2 /N 2 + 5v 3 + 10v 3 /N 2 

7 (1 - a) 7 + 21v(l - a) 5 + 70v 2 (l - a) 3 + 35u 2 (l - af/N 2 + 35w 3 (l - a) 
+70u 3 (l - a)/N 2 



Z L , a (N) = T ±- f n^e-^^^lrrnCl + ^-a), (13) 

where A L>a (N) = J ]jf =1 d^ i exp-f Tr ^ v ^ i ^. The general form of Z Li(X (N) for the matrix 
model with linear interaction from Eq. (12) using Wick theorem is 



Z L>a (N) = (l- a ) i + (l-a)( i - 2 )^^ + (l-a)( L - 4 ) £ VyVu 

i<j i<j<k<l 

+(l-a)( i - 4 ) Y, VuV jk +(±f(l-a)^ Wfi 

i<j<k<l i<j<k<l 

+.... (14) 

The partition functions for different a at a particular L can be obtained from Eq. (14). This 
relation has been checked explicitly for L upto 7 from Eq. (2). Each term in the general partition 
function Eq. (14) is accompanied by powers of (1 — a). The partition function corresponding to 
L = 4 (Table I) is given by Z L= ^ a (N) = (1 - a) 4 + 6u(l - a) 2 + 2t> 2 + v 2 /N 2 . In the Feynman 



diagram representation [ill ], this partition function represents a total of 10 diagrams (Fig. 2) where 
power of v gives the number of arcs, the coefficients of v give the number of diagrams with those 
many arcs and power of (1 — a) gives the number of unpaired vertices in the diagram. Such terms 
correspond to planar diagrams. The terms with powers of 1/N 2 represent non-planar diagrams 
(with crossing arcs) which have non-zero genus. The first term in Zl=4 i01 (N) is a planar term 
with no arcs (v°) and each unpaired vertex (4 in number) associated with a factor (1 — a), the 
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FIG. 2: The figure shows all possible diagrams for an RNA chain with length L = 4 in the linear interact- 
ing matrix model of RNA. The paired vertices are connected by dotted arcs and each unpaired vertex is 
associated with a factor (1 — a). 



second term represents 6 diagrams with one arc each (v 1 ) and each unpaired vertex (2 in number) 
accompanied by (1 — a), the third term corresponds to two diagrams with two arcs each (v 2 ) and 
no unpaired vertex at all and the last term represents one diagram with two crossing arcs (v 2 and 
1/N 2 ) i.e., a tertiary term with genus one and no unpaired vertices. For a = 1, the partition 
function for L = 4 is given by ZL = ^ a= \{N) = 2v 2 + v 2 /N 2 (Table I) i.e., a total of 3 Feynman 
diagrams (in Fig. 2, the diagrams corresponding to v 2 for g = and g = 1). For this a value, 
only those diagrams with completely paired vertices remain. Thus the structures in the model can 
be separated in two regimes as (i). < a < 1 which comprises of Feynman diagrams with both 
paired and unpaired vertices and (ii). a = 1 with only completely paired vertices. In (i). different 
genus structures with different weights (powers of (1 — a)) associated with each unpaired vertex 
((1 — a)) for different a's are included. The completely paired structures in (ii). consist of only 
those structures which have no unpaired vertices at all. 

If now = 1 is considered, the general partition function Eq. (14) becomes Zi t0l (N) = (1 — 
a) L + {1 _ a) L-2 Ei<j v .. + {1 _ a) L-A Y, i<j<k<l Vij V kl + (l-a) L - 4 E i<j<k< i V lk V jl + .... For an RNA 
chain in three dimensions, the partition function can be written as Z = f rifc=i ^ 3r fc/({ r })^L({ r }) 
where r k gives the position of the k-th. nucleotide in the chain, f({r}) is a model dependent function 
of geometry of the molecule and accounts for steric constraints of the chain, Zi{{r}) is the partial 
partition function where pairing between nucleotides comes with a factor Vij(rij) and the coupling 



constant of each unpaired and paired nucleotide is 1 [6|, 111]. If all steric constraints [/({^})] and 
spatial degrees of freedom (r^) are neglected and the chain is assumed to be infinitely flexible, 
the partition function becomes Z = Z L = 1 + Yli<j v ij + Y<i<j<k<i Vij v ki + Yli<j<k<i VikVji + •••• 
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FIG. 3: The figure plots genus distribution functions when a = 1 for the weighted normalized diagrams 
aL,g,a/Na with (a), length L keeping genus g fixed (g — 0, 1,2,3,4, 10,20), and (b) genus g with L fixed 
(L = 40, 80, 120,160). 

where Vij = e~^ e ' 1 '. The equation differs from the matrix model partition function Eq. (6) in 
[(J or Eq. (14) in the powers of 1/N 2 which classify terms with different topological character. 
Now, if it is assumed that each unpaired nucleotide in Zl is associated with a weight y, then 

Z = Z L = Z L>y = y L + y^ Y,i<j Vij + y (i " 4) Ei<j<k<l VijVkl + Ei< j<k <i V lk V 3 i + 

The external linear interaction is found to produce this weight on all unpaired bases in partition 
function of the chain, Eq. (14). Therefore, one can relate the partition function of the linear 
interacting random matrix model to that of an RNA with L nucleotides in three dimensions where 
each unpaired base comes with a weight y = (1 — a). 



B. Genus distribution functions 



The genus distribution functions are studied to observe the effect of external linear interaction 
on the distribution of pseudoknots among the structures of RNA. Figure 3 shows the genus dis- 
tributions for a = 1 where Fig. 3(a) plots the weighted normalized diagrams a,L,g, a =i/-^'a verses 
L at fixed genus and Fig. 3(b) shows the weighted normalized diagrams a,L t g ta =\/N a verses g at 
a particular L. In Fig. 4 the genus distributions are plotted for a = 0,0.75,1 for: (i). a pair 
of successive even and odd length L = (10, 11) by varying g [Fig. 4(a) and Fig. 4(b)] and (ii). 
different L at a fixed g = (0,3) [Fig. 4(c) and Fig. 4(d)]. For a chosen L, plots for a = 0,0.75,1 
are compared in Fig. 4(a) and Fig. 4(b). In the odd L = 11 plot, the a = 1 curve is absent (due 
to the absence of partition function for odd L's at a = 1, Table I). It is observed that the curve 
corresponding to a = 0.75 in the L = 11 plot comprises of points which are an average of points 
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FIG. 4: Genus distributions for the weighted normalized diagrams aL^, a /M a are plotted for a = 0,0.75, 1 
when (a). L = 10 and (b). L — 11 at a fixed genus and (c). g — and (d). g — 3 at a fixed length. For odd 
L (L = 11) the a = 1 distribution is absent. 

in the a = 0.75 and a = 1 curves when L = 10 (same genus points for the two L are considered). 
For example in Fig. 4(a), the points corresponding to g = 2 for a = 0.75 and a = 1 curves 
are averaged and the cll q a /M a value obtained is similar to the aL, g , a /M a value for g = 2 in the 
a = 0.75 curve in Fig. 4(b). The same is observed for other such even and odd length comparisons 
and also when aL, g , a /M a is plotted with L for fixed genus, g = (0,3) (Fig. 4(c) and Fig. 4(d)). 
Notice that for a = 0.75 in the fixed g plots, the points for each successive even and odd L lie 
close together separated by large distances from the nearest neighboring even and odd L pair. The 
values of the weighted normalized distribution can be found from these figures, for example, in 
Fig. 4(c) for a L=4i9=0 ,a=o = 9 and M a=0 = 10, 0^=4,9=0,0=0/^=0 = 0.9 and for a L=4iS= o, a =i = 2 
and M a= \ = 3, a_L=4 iS =o,a=i/A/" a =i = 0.66. It is also found numerically that the maximum value of 
genus g for the matrix model with linear external interaction is g < L/4 for even L and all a and 
for odd L and < a < 1, while for a = 1 there are no structures. 

A small length even and odd behavior of the partition function starts emerging as a is increased 
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from towards 1 and becomes very distinct at a = 0.75. This behavior is more and more prominent 
for a's greater than 0.75 and in the close vicinity of a = 1. The random matrix model with a linear 
external interaction therefore shows small differences in the even and odd L partition function at 
small L. The length scale in the linear interacting matrix model of RNA, below which the even 
and odd length partition functions are distinguishable, depends upon the interaction parameter a 
and is found to be L a ~ 1/(1 — a) 2 . 

III. SOME INTERESTING RESULTS IN THE RANDOM MATRIX MODELS OF RNA 

A. Scaling relation between the linear interacting and re-scaled random matrix models of 

RNA 

Starting with the generalized partition function Zi ja (N) equation (14) and factoring out (l—a) L 
from each term results in 

Z L (v,a,N) = (l_ a ) L [l + (l-a)- 2 ^K i + (l-a)- 4 £ V i5 V u 

i<j i<j<k<l 

+(l-a)~ 4 £ V u V, h + (±) 2 (l-a)-* Y, ^ikVji 

i<j<k<l ' i<j<k<l 

+....]. (15) 

The partition function is written with v and a in the argument as its dependence on these 
parameters will be discussed here. The bracketed quantity on the right hand side of Eq. (15) is 
the generalized partition function of the random matrix model of RNA (a = 0) where base pairing 
interaction strength v has been re-scaled by (1 — a) -2 . The general partition function of the matrix 
model of RNA with a linear external interaction can thus be written as 



Z L (v,a,N) = (l-a) L Z L [- 



r AN]. 



(16) 



(1-a) 2 ' 

This is the scaling theory for the partition function of the random matrix model of RNA with 
linear external interaction with scale factor (1 — a) L having a scaling exponent L and a crossover 



exponent^ Such sea 
systems 



21 



23, 



ing relations have been found in other systems like the anisotropic magnetic 



24j (section 4 of the first reference in [22], equations (6-10) of the second 
reference in 22] and equation (5.4) in 23]) but in the context of random matrix models such a 
scaling form has not been observed before in the literature (to the best of our knowledge). For 
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the a = partition function, which corresponds to the a = line in Fig. 1(a), re-scaling v by 
(1 — a)~ 2 and multiplying by (1 — a) L will give the other q / lines in Fig. 1(a) as well as the 
a = 1 line where Zl a(N) = for odd lengths. The distribution functions (Fig. 3 and Fig. 4) can 
be equivalently studied from the re-scaled random matrix model of RNA for different a values. The 
matrix model of RNA with linear external interaction can thus be viewed in two equivalent ways, 
(i) . the linear interacting matrix model of RNA where addition of a linear external interaction in 
the action of the partition function [11] results in each free base of the chain getting weighed by 
(1 — a) and (ii). the random matrix model of RNA where the base pairing strength v is re-scaled by 
a factor (1 — a)~ 2 with an overall scale factor (1 — a) L . In the present work, the first view point is 



followed as it incorporates t 
formalism presented here 



le study of many more general interactions and their effects using the 



33J. The addition of more complicated interactions in the matrix model 



action in general, will not produce such a scaling relationship. 



B. Thermodynamics 

Thermodynamic properties such as the free energy and specific heat for the linear interacting 
matrix model or equivalently the re-scaled random matrix model of RNA are calculated as functions 
of length of the chain L (which may be useful in stretching experiments), a and temperature T 
from Eq. (12). The partition function Zl iCi (N) depends upon temperature T through v which is 
given by v = e _/3e . The Boltzmann constant, fee, and e are set equal to unity. The calculations 
are performed with T = 2 (except for where the T dependence is studied) and N = 100 and it 
is expected that the analysis will go through for very large N. In the re-scaled random matrix 
model of RNA, the a = 1 line corresponds to the limit v' — ► oo as v' is given by v' = v/(l — a) 2 



Thermodynamic properties (free energ y, c hemical potential and specific heat) of the model in 
have been calculated and discussed in 10| as a function of temperature T. 



11] 



1. Free energy and phase transition 

The free energy F is found numerically from the partition function ZL a (N) Eq. (12) using the 
relation Fl i0i (T,N) = — ^ln[ZL ja {T, N)]. As a function of T, the free energy is shown in Fig. 5 
for L = 19 and L = 20. The curves for different a, for both even and odd lengths, lie above the 
preceding smaller a as a — ► 1. The free energy curve for a = starts from zero at low temperatures 
and remains negative as T is increased. For < a < 1, F at small T rises from zero to a peak +F 
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FIG. 5: The figure plots free energy as a function of temperature T for (a). L — 19 and (b). L = 20 (the 
a = 1 curve is absent in (a) as F is imaginary). Plots in (c,e) show F with L and a in a 3D plot. Plot 
(c) shows all lengths (both even and odd) where the small L even and odd behavior is seen, plot (d) shows 
only the even lengths and plot (e) shows F for odd lengths only. The F verses a plot in (f) for L = 19 
shows a cusp (inset) as a — ► l - . The symbols in plots (a) and (b) correspond to the following a values: Box 
— > a = 0, Triangle — > a — 0.25, Star — > a — 0.5, Diamond — > a = 0.75 and Pentagon — » a = 1. 

value and falls down before crossing over to the —F region. The a = 1 curve (in the L = 20 plot) 
starts from a non-zero value and drops gradually to the — F side. The a = 1 curve is absent in 
the L = 19 plot as i* 1 becomes imaginary. The positive to negative crossover of F for different a is 
at different T given by T a ~ a. This identifies a low and high temperature behavior of F with T 
for different a's. Three dimensional plots of F with length of the chain L and a are shown in Fig. 
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FIG. 6: The figure shows first derivative of free energy, (dF/da), with respect to a for a pair of odd 
(L = 9, 19) and even lengths (L = 10, 20). The numbers in bracket of the legend of each plot correspond to 
the numerical a = 1 values. 

5(c)-Fig. 5(e). The small L even and odd behavior of the partition function (discussed in Section 
2.2) is visible in Fig. 5(c). Figure 5(d) shows F plotted with only the even lengths and a. At small 
lengths for a = 1, the plot is flatter compared to relatively large lengths where it is sharply peaked. 
The peak value however becomes more negative with increasing L. It is expected that as length 
is increased further, the peak will become more and more sharp near a = 1. This gives signature 
of a phase transition, at very large L, at the a = 1 line for the even lengths. The free energy for 
odd lengths and different a in Fig. 5(e) shows that F becomes imaginary at a = 1. The peaked 
behavior is present for all, small and large, L with the peak getting sharper for larger lengths. In 
Fig. 5(f), the cusp from the a — > 1~ side for F is displayed for a particular L = 19. Hence, for odd 
lengths the free energy exhibits a phase transition in the linear matrix model of RNA. In order to 
understand the transition better, the first derivative of free energy with respect to a as a function 
of q is shown for different L in Fig. 6 (Fig. 6(a) for L = 9, 19 and Fig. 6(b) for L = 10, 20). For 
the odd length plot in (a), the L = 9 and L = 19 curves rise smoothly as a is increased from 
towards 1. But in the very close vicinity of a = 1, the curves show a steep rise with the values at 
a = 1 being very large (Fig. 6(a)). The curves however are seen to pass through zero very close to 
a = 1. In the even length plot (b), the L = 10 and L = 20 curves behave similar to plot (a) for a 
from to near 1 but suddenly goes down to approximately zero very close to a = 1. These plots 
indicate a different behavior for the odd and even lengths near a = 1. The scaling found in Eq. 
(18) and discussed in Fig. 1(a) is also present close to the transition line a = 1. 
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FIG. 7: The figure shows specific heat C v , for a = 0, 0.25, 0.5, 0.75, 1, plotted as a function of (i). temperature 
T for (a). L = 19 and (b). L = 20 (with a = 1 plot in the inset) and (ii). length L. The different a curves 
lie one above (below) the other when plotted as a function of T (length) as a — > 1. The symbols in the plots 
correspond to the following a: Box — > a = 0, Circle — > a = 0.25, Triangle — > a = 0.5, Star — > a = 0.75 and 
Pentagon — > a = 1. 

2. Specific heat 

The specific heat can be found from the free energy as C v = —Td 2 F/dT 2 where constant volume 
implies constant L here. The C v verses T plots for different q(= 0, 0.25, 0.5, 0.75, 1) and a pair of 
odd (L = 19) and even (L = 20) lengths is shown in Fig. 7(a) and Fig. 7(b) respectively. The 
curves for different a lie above the preceding a curve as a — > 1. The maximum C v for different 
a curves increases as a is increased towards 1. For L = 19, the a = 1 curve is absent as free 
energy is imaginary (as discussed above). The a = 1 plot for L = 20, shown as an inset to Fig. 
7(b), is indeterminate at T = and oscillates wildly at low T. The scale of the inset (a = 1) 
is of the order 10~ 12 smaller than the other a / 1 curves. As a function of L, C v is plotted for 
different a values in Fig. 7(c). At small lengths, the even and odd behavior (as found in the genus 
distribution functions) can be observed in this plot as a approaches 1. The different a^0 curves 
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FIG. 8: The first derivative of specific heat, (dC v /da), with respect to a for a pair of odd (L = 9, 19) and 
even lengths (L — 10, 20) is shown. The numbers in bracket of the legend (in each plot) correspond to the 
numerical a = 1 values for the particular length curves. 

lie below the preceding a curve. In the a = 1 curve, C v for odd lengths is indeterminate and is 
zero for even lengths. To explore the behavior near a = 1 line, the first derivative of specific heat 
with respect to a, dC v /da, plotted as a function of a is shown in Fig. 8 (Fig. 8(a) for L = 9, 19 
and Fig. 8(b) for L = 10,20). In the odd length plot (a), the L = 9 curve bells down very near 
a = 1 and the L = 19 curve shows a rising behavior with very large values in the close vicinity of 
a = 1. The curves however pass through zero very near a = 1. The plot (d) shows a smooth rise 
initially which shoots upward sharply near a = 1 and drops to approximately zero at a = 1. This 
again indicates a different behavior for the odd and even lengths near a = 1. 

Therefore, it is observed that as the length is increased, the plots (for both F and C v ) become 
more sharply peaked near a = 1. It is expected that as one goes to very large lengths (of the order 
of hundreds of bases), the derivatives with respect to a as a function of a will be discontinuous as 
they approach (i). +oo and — oo for dF/da and (ii). — oo and +oo for dC v /da when a — > 1~ and 
1 + respectively (for the odd lengths see also Fig. 5(f)). The transition for the odd and even lengths 
seem to be different. In order to completely establish the phase transition in these matrix models 
of RNA and to make the numerical results presented here more concrete, a detailed analytical 
study for small L and near a = 1 is required. Moreover, the appearance of a small L even and odd 
behavior in the distribution functions, free energy and specific heat needs further understanding 
since these matrix models are models of RNA and hence this phenomena is a possible prediction 
for real very small RNAs (very little experimental work has been done in this direction). 
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IV. CONCLUSIONS 

An analytic calculation of the partition function after introducing a linear external perturba- 
tion in the action of the partition function of the random matrix model of RNA in [ll| is given. 
A perturbation parameter a is found which can be tuned to different values to get two different 
structural regimes in the model. The presence of different regimes was seen in the study of (sec- 
ondary structures of the) RNA-like polymer models 25[]. In this work, different regimes are found 
with tertiary structures taken into account in addition to the secondary structures. The structural 
changes between two regimes with respect to a can be defined as comprising of (i). unpaired and 
paired base structures, in the regime < a < 1 and (ii). completely paired base structures when 
a = 1. The regime (i) is like the molten phase of the RNA structures shown by many statistical 



models of RNA (see for example [261]). The fully paired regime corresponds to a range of struc 



tures (similar to the fully paired regime discussed in 25|). The general form of Z^^N) has terms 
weighted by powers of (1 — a). The linear perturbation Wi = w introduced in the action of the 
partition function in Eq. (2) thus alters the contribution of unpaired bases as compared to the 
paired ones. In the linear interacting model picture, a physical effect which gives rise to w can 
be the different surroundings of the cell in which the RNA lives such as the presence of ions or 
molecules 27J, |28I . |29|, then (1 — a) represents this change. The interaction introduced here may 
have other concrete realizations, for example in the pulling experiments ([f], [3^] and references 
therein). The genus distributions at different L, a and genus g exhibit small differences for, even 
and odd, short lengths with these differences disappearing at long lengths in the partition function. 

It is observed that the linear interacting matrix model is related via a scaling expression to the 
random matrix model of RNA where the base pairing interaction strength v has been re-scaled by 
(1 — a)~ 2 . The two approaches are equivalent in this sense. The thermodynamic variables (free 
energy and specific heat) as a function of (i). L show a different even and odd length character 
present in the model, which may be of interest in the stretching experiments, (ii). T show a 
different low and high temperature behavior in the free energy. The analysis of the free energy of 
the interacting matrix model with respect to a shows a phase transition for the odd lengths and 
a sharp transition is expected for the large even lengths. A detailed analytic study of the phase 
transition reported here (found numerically) in these matrix models of RNA is for a future work. 

The formalism gives a systematic analytical method of studying the effect of external interac- 
tions on the number of the planar and non-planar diagrams for a given L and on their thermody- 
namic properties. These models are matrix models of RNA which enumerate structures of all types 
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that are possible for a given length of the chain. So the characteristics observed in the distribution 
functions and the thermodynamics, like the small L even and odd behavior, should be found in 
some real RNA's such as the micro-RNA which are small with only 21-23 nucleotides. The results 
found here are also interesting in their own right purely from the point of view of random matrix 
theory. 
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